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ABSTRACT 

So far, six mechanisms have been proposed to account for the Galactic disk heat- 
ing. Of these, the most important appear to be a combination of scattering of stars 
by molecular clouds and by spiral arms. We study a further mechanism, namely, the 
repeated disk impact of the original Galactic Globular Cluster population up to the 
present. We find that Globular Clusters could have contributed at most a small frac- 
tion of the current vertical energy of the disk, as they could heat the whole disk to a z 
— 5.5kms _1 (c.f. the observed 18 and 39 kms -1 for the thick and thin disks respec- 
tively). We find that the rate of rise of disk heat (a=0.22 in a z ~ t a with t being 
time), is close to that found for scattering by molecular clouds. 

Key words: galaxy: disk - globular clusters: general. 



1 INTRODUCTION 

The disk stars in the Galaxy can be described by a DF rep- 
resented by a Schwarzschild distribution where the number 
of stars with velocity v within a small range d z v is: 
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where no is the number of stars per unit volume, and the 
(tr , cr^ , a z are the dispersions in the three cylindrical coor- 
dinates (Binney & Tremaine 2008). This ellipsoid's shape is 
similar for all stars in the disk, but its size depends on the 
stellar population. For example, the observed dispersion is 
three times larger for cool, red stars than for hot blue stars. 
This suggests that stars are born on nearly circular, cold, or- 
bits that heat up as time goes by. The value determined for 
the current vertical heating of the thin disk is 18kms _1 , that 
for the thick disk is 39kms _1 (Binney & Merrifield, 1998). 

The nature of the heating agent is discussed in detail 
in Binney & Tremaine (2008), who list the following pos- 
sible candidates: encounters with MACHOs from the dark 
halo, velocity kicks from passing stars, scattering by spiral 
arms, scattering by molecular clouds, mergers with satellite 
galaxies, and finally, substructure in the dark halo. They 
conclude that the most likely mechanism is a combination 
of spiral transients and scattering by molecular clouds. 
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In this paper we examine whether the repeated impact 
of Globular Clusters (GC) since Galaxy formation could be 
a significant mechanism for disk heating. In Section 2 we 
examine the physical processes involved. Section 3 presents 
numerical results for the current GCs. Section 4 shows the 
results for the entire GC population since Galaxy formation, 
(1) then Section 5 discusses these results. 



2 PHYSICAL PROCESSES 

The Galaxy formed ~ 10 10 years ago, with its population of 
Globular Clusters moving along orbits in the Galactic poten- 
tial. Assuming that the initial disk of the Galaxy was cold, 
the GC impacts through the disk will gradually heat up the 
disk as a result of dynamical friction. This is a process where 
a heavier object of mass M gravitationally interacts with a 
large number of randomly distributed less massive objects 
of mass m, and transfers some of its kinetic energy to the 
less massive particles. The expression for the deceleration 
was first given by Chandrasekhar (1943), and re-expressed 
in a convenient form by Binney & Tremaine (2008) as: 
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where vm is the velocity of a GC of mass M, the individ- 
ual disk stars having a massm <g M, In A is the Coulomb 
logarithm, whose value is set to 18.1, the value in the Solar 
neighbourhood (Binney & Tremaine, 2008). The parameter 
X = vm I (\/2<r), and a is the velocity dispersion of disk 
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stars. The local mass density is p. Equation 2 assumes that 
the heavier object is a point mass, but tests have shown 
that it is nevertheless a remarkably good description of the 
drag experienced by a finite body orbiting in a stellar sys- 
tem (White, f978). As a result, the GC experiences a force 
directed opposite to the direction of motion, along a path 
of length A of disk material equal to the GC pathlength 
through the disk. The energy transfer Et is then: 
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so that the GC is slowed down and moves to a lower orbit. 
Ultimately, the orbit will decay to the point where the GC 
dissolves into the bulge or disk. Formally, Equations 2 and 
3 are obtained for an infinite uniform medium, so represent 
an approximation for a GC crossing the Galactic disk. 

The dynamical friction experienced by a GC in the halo 
is small compared to that during its transit through the disk. 
For example, in the case of uiCen, the dynamical friction 
from a disk transit is two orders of magnitude greater than 
that it experiences during the halo part of its orbit. 

With time, several phenomena occur, which strip the 
GC of mass: stellar evolution which drives mass out of the 
GC via winds and supernovae, two-body relaxation where 
the cumulative effect of gravitational interactions between 
stars occasionally provides sufficient energy for a star to es- 
cape the cluster, and tidal shocks, where as the GC tran- 
sits through the disk, it experiences a rapidly changing tidal 
force which can increase the stars' kinetic energy to the point 
that they escape the cluster. The result of these phenomena 
on the Galactic GC population has been studied extensively 
by Fall & Zhang (2001), and Prieto & Gncdin (2008). The 
latter study the dissolution of GCs in a constant combined 
halo-disk potential for GC masses in excess of 10 5 M Q with- 
out orbital decay. The former work back from today's Galac- 
tic GC population, to predict the GC-IMF, based on a static 
spherical potential, and accounting for GC mass loss. 



3 HEATING DUE TO THE EXISTING 
GLOBULAR CLUSTERS 

There are ~160 GCs in the Galaxy today, with 
most listed in the Harris catalogue (Harris, 1996), 
of which the latest on-line version is available at 
www.physics.mcmaster.ca/Globular.html, dated 2003. In 
this section we describe how we calculate the disk heating 
were it caused exclusively by these 160 GCs. To do such a 
calculation requires a knowledge of the GC trajectories, so 
that energy transfers at disk crossings may be evaluated. 
Determining trajectories in turn requires a knowledge of 
GC positions and proper motions. Of the 160 GCs, only a 
subset of fifty-four have reliably determined proper motions 
(Moreno & Pichardo 2006 and 2008, and Kalirai 2007). We 
start by determining the heating due to this subset of fifty- 
four GCs. It is possible to determine their past trajectories 
up to the formation of the Galaxy ~ 10 10 years ago if we 
assume that the Galactic potential has remained constant 
throughout the time since Galaxy formation. Note that the 
disk potential is assumed constant, with the dissolving GCs 
having no influence on disk mass. We can then establish the 



characteristics of their crossings over that time (position in 
disk, velocity, time, angle of incidence on disk). 

We follow for the orbit calculations the method de- 
scribed in Vande Putte & Cropper (2009), using the poten- 
tial proposed by Fellhauer (2006). This potential has three 
components, representing a bulge, halo, and disk: 



(4) 



The contribution of the bulge is given by a Hernquist 
(1990) potential: 
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where M b = 3.4 x 10 10 M Q , and a = 0.7 kpc. 

The axisymmetric disk potential is represented by a 
Miyamoto-Nagai (1975) expression: 
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with M d = 10 11 M Q , a = 6.5kpc, b = 0.26kpc. A 

spherically-symmetric halo contributes to the potential in 
a logarithmic form: 

*h = f HR 2 +* 2 + d 2 ) (7) 

with vo — 186kms _1 , and d = 12kpc (Fellhauer et al, 2006). 

In our orbit calculations we track the GCs, and note the 
crossings data such as time, position of crossing, and angle 
of incidence on the disk. In doing this, we found that the 
fifty-four GCs cross the disk plane 10454 times over the last 
10 Gyr. We assume an initial disk star velocity dispersion in 
the z direction, of lkms -1 , based on the study of vertical ve- 
locity dispersion in disks as a function of time, by Jenkins & 
Binney (1990), thus finding the initial vertical energy in the 
disk. We then use Equation 3 to calculate the energy transfer 
to the disk at each crossing. Vande Putte & Cropper (2008) 
show how the track of the GC through the disk focusses 
stars towards it, essentially randomising the velocities, so 
we divide this energy by 3 and add the already-present ver- 
tical energy to find the total vertical heating of the disk. We 
assume instantaneous spreading of the energy throughout 
the star population in the disk. When applying Equation 3, 
we combine p and A to give the surface density E and we 
consider the following relation for the surface density of the 
stellar disk (Binney & Merrifield, 1998): 



E(i?) = So exp(-R/R d ) 
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where R is the radial distance of the impact point, and 
Rd = 2.5kpc (Binney & Tremaine, 2008). Combining with 
the requirement for a total disk mass of 4.5 x 10 10 Mq (Bin- 
ney & Tremaine 2008) gives: 



p(R) = 4.46exp(-i?/2.5kpc) M Q pc~ 



(9) 



Widrow, Pym & Dubinski (2008) find independently 
thati? d = 2.8± , ; 2 Jkpc and M d = 4.1±g;|| x 10 10 M Q . These 
values are similar to those used in the potential above. 

As the visible disk is about 15kpc in radius, the calcu- 
lation discards any crossing beyond that radius (there are 
only just over 200 such cases out of a total number of 10454 
crossings). 

To sumarise, the calculation proceeds with the following 
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steps. First, we calculate the disk's initial vertical kinetic 
energy based on the mass of stars and a velocity dispersion of 
lkms -1 (this is 4.5 x 10 j3 ergs) . We then calculate the energy 
transfer Et for the first crossing, using Equation 3, and 
spread the transferred energy uniformly over all stars in the 
disk and calculate a z , the vertical velocity dispersion based 
on 1/3 of the energy transfer being available for vertical 
heating. Finally, we repeat the steps above for all crossings 
in turn, for all fifty-four GCs. 

This leads to a final, present-day a z — 1.8kms _1 from 
all the crossings by the subset of fifty-four GCs. Running the 
calculation with wCen only, the most massive GC, produces 
a z = 1.6kms _1 , which indicates the importance of GC mass 
in the energy transfer. Other factors, related to the GC or- 
bit, influence the energy deposition, mainly the number of 
crossings, and their location relative to the Galactic cen- 
tre. For example, 47Tuc is the second most massive GC in 
our subset of 54, being a factor 3.3 lighter than wCen. On 
mass alone, Eqn.3 predicts an energy deposition for 47Tuc, 
a factor about 11 less than for uCen. In fact, we calculate 
the energy deposition to be a factor 75 less. This is because 
47Tuc has only 138 crossings, whereas wCen experiences 262 
crossings. In addition, the radial position of uCen crossings 
ranges from lkpc to about 6kpc, with a quarter being near 
lkpc. On the other hand, 47Tuc crossings range from just 
under 6kpc to just over 7kpc, i.e. in a less dense region. 

To obtain the heating from all 160 GCs, we make the 
assumption that the 54 orbits and masses are representative 
of all 160 orbits and masses, and multiply the energy transfer 
at each crossing in the above calculation by the ratio of 
160 to 54. The result is a present-day disk heating of o z — 
2.8kms~ 1 . This is clearly an overestimate, as it implies that 
there are ~ 3 wCen equivalents among the current total of 
160 GCs. Nevertheless, this result shows the relatively small 
effect the current population could have had on disk heating. 

We have examined the consequence of regarding the GC 
as a point mass, whereas in reality it is a system with many 
degrees of freedom. We follow Binney & Tremaine (2008) 
and find that the total energy transferred by disk crossing 
to the stars in the GC is 

AE = (1Q) 

where E d is the surface density of the stellar disk, and a 
the average semi-major axis of stars within the GC. The 
quantity in square brackets in Eqn.3 is close to unity in 
all cases considered here. Hence, the ratio of the dynamical 
friction (Eqn. 3) to the quantity represented by Eqn.10 is 
then ~ > as p A = E^. With a disk surface density 

of 45M pc~ 2 , a ~ 1 pc and the mass of those GCs that 
play the overwhelming part in energy transfer > 1O 5 M , the 
ratio of the energy transfer to the disk by dynamical friction 
compared to the energy transferred to the stars within it is 
>10 4 . Hence the fraction of the orbital energy diverted into 
tidal effects is negligible. 

4 EFFECT OF ALL GCS SINCE GALAXY 
FORMATION 

While the fraction of the orbital energy diverted to tidal ef- 
fects is a small fraction of the orbital energy, it is a larger 



fraction of the binding energy. Hence the disk transits and 
other effects in Section 2 contribute to the dissolution of 
GCs. The number of GCs in the past was therefore much 
higher and we now consider the effect of this on GC heating 
contribution. We follow all GCs ever present from ~10 10 
years ago, to the present, using mass functions from Fall 
& Zhang (2001). We select a Schechter function from their 
figure 3 for the GC-IMF. This has the form xjj {M) oc 
M 13 exp(-M/M«) with (3 = -2 and M„ = 5 x 10 6 M . 
We chose this IMF because it yields a close fit to the ob- 
served distribution at the current epoch, while including a 
large component of initially present GCs. Specifically, using 
the data in Fall & Zhang, the GC-IMF is: 

MM) ^ 2 x 10 8 M~ 2 exp(-M/M„) (11) 

and the evolved mass functions are given for t=1.5, 3, 6, 
and 12 Gyr. These have a more or less lognormal form, with 
peaks gradually moving to higher masses. 

We again make the assumption that the orbits and 
crossings of the fifty-four GCs are representative of all GCs 
ever present. It is then possible to calculate the disk heating 
for, say, N GCs with this mass spectrum, by assuming that 
for each crossing by one of the fifty-four, the energy transfer 
is in fact N/54 times larger. 

In the calculations we use 16 bins to cover the mass 
range 10 3 to 1O 7 M . The total initial number of GCs is 
~160,000. We account for GC mass loss, hence varying mass 
function, by applying the evolving mass functions of Fall & 
Zhang (2001). As these authors use only four time subdi- 
visions, we interpolate between their curves, for each mass 
bin, to obtain 14 time subdivisions. We calculate the heat- 
ing due to the complete history of heating caused by these 
populations depositing energy in the disk through dynam- 
ical friction for each energy bin and time interval before 
summing to obtain the total energy input. 

The above calculations lead to our finding that GCs 
could on their own heat the disk to a value a z — 5.5kms _1 . 
This value lies below the thin and thick disk values of 18 
and 39kms~ 1 respectively. 

These calculations intrinsically carry the time history of 
the a z dispersion. We have therefore plotted the dispersion 
as a function of age in Figure 1. Jenkins & Binney (1990) 
found that the dispersion increases with age as a z ~ t a . 
We find from fits to our data that a=0.22. We are able 
to compare this with observational data such as that from 
Nordstrom et al (2004) who performed a kinematic analysis 
of over 16,000 nearby F and G stars, and who found a value 
of a~ 0.47. 

To obtain a value of a z — 26kms _1 , a value predicted 
for a population aged lOGyr from the Nordstrom data, the 
energy input from GCs would need be twenty five times 
larger than provided here by the GCs. In other words, the 
GCs are at most a minor contributor to disk heating. 

Our value of a z can be seen as an upper estimate, be- 
cause the radial extent of the disk was smaller in the past, 
and because the number of stars was smaller in the past, two 
facts unaccounted for here. On the other hand, orbital de- 
cay due to dynamical friction has been neglected, thus also 
overestimating the energy transfer calculated here, but this 
neglect has also underestimated the initial GC-IMF. 

We have checked the possible contribution of stars 
stripped from the GCs by the effects in Section 2 to the 
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Figure 1. Velocity dispersions ay/ for single stars with relative 
age errors < 25% as a function of age from Nordstrom et al (2004, 
solid circles), with our corresponding values of <r z due to GC 
impacts, cry and at the thick and thin disk values. 



vertical velocity dispersion of the disk. The total mass of 
stars from the evaporated GCs using Fall & Zhang's GC- 
IMF is ~ 1.5 x 10 9 M Q , or <4% of the total disk mass. 
Because the fraction residing in the disk at any instant is at 
least an order of magnitude smaller, even if these stars are 
kinematically hot, their contribution can be neglected. 

The particular selection of GC mass functions may have 
a bearing on the result. Parmentier & Gilmore (2007), for 
example, discuss bell-shaped mass functions. These are par- 
ticularly at variance with the GC-IMF we selected from Fall 
& Zhang, which rises monotonically as the mass decreases. 
To investigate the impact of this, we assume that the GC- 
IMF we chose, drops precipitously to zero near the peak of 
the current mass distribution, towards lower values of mass. 
We calculate that the energy input would drop only 4%, 
with almost no effect on a z . This underlines the importance 
of the massive clusters. 



disks (18 and 39 kms -1 , respectively). Using the vertical 
velocity dispersion of Nordstrom et al (2004), the GC ver- 
tical heating contribution is about one twenty fifth that of 
the total disk vertical heating, making GCs a second order 
contributor to Galactic disk heating. As the energy transfer 
scales with mass squared, an extension of this latter result 
has the implication that the upper limit to the total mass 
associated with a GC is about five times the stellar mass. 
This places a constraint on the possible dark matter content 
in GCs, as far back as the time of Galaxy formation. 
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5 CONCLUSION 

We have investigated the vertical heating of the Galactic 
disk's stars, due to GC crossings inducing dynamical fric- 
tion. We achieved this by modeling the orbits of a subset 
of fifty-four GCs with well-determined orbital parameters, 
and by assuming a static Galactic potential and disk stel- 
lar population. Assuming these orbits and masses are rep- 
resentative of all GCs at all times, these results were first 
extrapolated to the 160 existing GCs to find their effect on 
disk heating. Next, to consider the effect of all GCs ever 
present, we used time-dependent mass functions represent- 
ing all GCs ever present in the Galaxy. The resulting value 
of o z = 5.5kms _1 is smaller than that of the thin or thick 



